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1. Introduction 

Scaling concept plays a vital role in describing critical phenomena associated with a 
second-order phase transition. The fundamental hypothesis states that, in the vicinity 
of a critical point, the largest length scale of the fluctuation of the order parameter 
diverges and all length scales contribute with equal importance E]. Theoretical 
arguments based on this hypothesis explain why most thermodynamic quantities near 
the critical point exhibit power-law behaviour with characteristic exponents that are 
independent of the microscopic details of a system jSllllEl- 

A primary example of a physical model exhibiting a second-order phase transition is 
the two-dimensional Ising model with ferromagnetic interaction [HI 13 IE]- This model has 
been used extensively for innumerable projects in statistical physics, mainly due to its 
simplicity and broad applicabihty to real systems. Whereas its physical properties have 
been thoroughly investigated, it still continues to raise interesting issues that are relevant 
to a wide range of critical phenomena. Most intriguing among them is the critical 
properties of the Ising model defined on curved geometry UHl IHl 1121 HH Ull Hi] • In 
fact, several studies have been carried out on Ising lattice models with topologies ranging 
from a torus [HI El HE] , and sphere jTU [TTl [T21 CHI Ull to genus two curved surfaces [12] , 
as well as those with Brascamp-Kunz boundary conditions fTlEl- The results of these 
studies are consistent with the fact that an alteration in the topology or the boundary 
condition of the Ising lattice does not change its scaling behaviour, and thus the system 
remains in the flat-space Onsager universality class 

The objective of the current study is to focus on an alternative property of curved 
geometries — surface curvature — instead of their topology. Our main concern is whether 
a uniform change in the constant surface curvature affects the scaling behaviour of 
the mounted Ising lattice model. It should be noted that in most systems considered 
thus far, the magnitude of surface curvature is spatially concentrated on a portion of 
the surface, even producing a conical singularity. This inhomogeneous distribution of 
local curvature would possibly make it difficult to distinguish the effect of curvature 
from among other incidental contributions on the critical properties. In addition, when 
a surface has a closed form (e.g., a sphere), its ability to attain the thermodynamic 
limit with a constant surface curvature is disabled |. This limitation can be removed 
successfully by employing a surface with a constant negative curvature. This surface, 
in which the Gaussian curvature possesses a finite constant value at arbitrary points, is 
simply connected and infinite [T^ 120] . Hence, such a surface can serve as an example 
for considering the geometric effects on the critical properties of the mounted system. 

In the present paper, we investigate the critical behaviour of the two-dimensional 
Ising lattice model defined on a curved surface with a constant negative curvature. 
Monte Carlo (MC) simulations and finite-size scaling analyses are employed to compute 
the critical exponent 7 for the zero-field magnetic susceptibility and fi for the correlation 

I For example, a sphere reduces to a flat plane at some point within the thermodynamic limit where 
the effect of curvature is absent. 
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volume. We demonstrate that the values of both 7 and /i deviate from those for the 
planar Ising lattice model, which indicates the relevance of the intrinsic geometry of 
the underlying surface to the critical properties of the mounted Ising model. Moreover, 
when reducing the boundary contributions, the values of 7 and /i exhibit a tendency to 
shift to those derived from the mean-field approximation. This non-trivial behavior of 
the critical exponents is qualitatively consistent with the conclusion based on the series 
expansion analyses and that deduced from the quantum field theory [22] ■ We also 
calculate the fourth-order Binder's cumulant that provides a check of the results of finite 
size scaling. 

2. Scaling arguments for the Ising model 

Let us briefly review the framework of the scaling argument that successfully explains 
the critical properties of the two-dimensional planar Ising lattice model jSj. The scaling 
hypothesis states that, in the vicinity of the critical temperature T^, the singular part 
of the free energy /s of the Ising lattice per site should be a homogeneous function: 

A'^/s(t,/i) = /s(A%A^/i), (1) 

where t = {T — T^) /T^, h = H/ {kBT) and H represents an external magnetic field. The 
parameter A'' with a spatial dimension d indicates the rescaling of the total number of 
sites from N to A'^A^; this results from the transformation of the linear dimension of the 
entire Ising lattice: L —>■ XL. By eliminating A from (^, we obtain the scaling relation 
of /s: 

/s(^^^) = l^r-^(^)> (2) 

where is a universal scaling function. The appropriate differentiation of Q yields 
the power-law form of thermodynamic quantities such as the zero-field susceptibility 
X = d'^fs/dh^ oc It]""^, where the critical exponents are expressed as functions of x and 
y (for instance, 7 = {2y — d)/x). 

Equation (PJ) is justified when the Ising lattice is defined on a fiat plane, since the 
rescaling of L by A is equivalent to that of by A"^. However, this is not the case when 
the Ising lattice is defined on a curved surface; while a wide range of regular lattices 
can be constructed on curved surfaces with a constant Gaussian curvature [1111201, the 
relation N = L'^ become invalid for these lattices due to the differences in the metric 
of the underlying geometry. Hence, the rescaling L XL does not imply — > A'^A^, 
which requires some modifications of (0). § 

We thus introduce an alternative rescaling parameter A for considering the scaling 
relation of the Ising model on curved surfaces: 

Afsit,h) = f,iA%Ayh), (3) 

§ Similar argument has been made regarding an infinitely coordinated Ising model; see Ref. [2^1 • 
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Figure 1. Schematic illustration of a regular heptagonal lattice in terms of the 
Poincare disk representation. The number of concentric layers of heptagons is r = 3 
in this figure. All heptagons depicted within the circle are congruent with respect to 
the metric given in Eq. ljA.3|) . The circumference corresponds to an infinite distance 
from the center of the circle. 

where A represents the rescahng of total sites AN. When the underlying geometry 
of the lattice is flat, the relation Q reduces to ((T)) since A = A'^. In this case, the 
parameters x and y given in (jH)) are defined as x = x/d and y = y/d, thus yielding 
identical values of critical exponents (for instance, 7 = {2y — l)/x = {2y — d)/x = j). 
In contrast, when the underlying geometry is curved, A can be no longer expressed 
as a power of A; thus, x and y are not related to x and y. Consequently, the critical 
exponents for the latter model, which are determined by x and y, may quantitatively 
differ from those for the planar Ising lattice. This naturally motivates us to evaluate the 
critical exponent directly by constructing the Ising lattice model on curved surfaces. 

3. Regular tessellation of curved surfaces 

A simple spherical surface seems to be the optimal geometry to consider the curvature 
effect on the critical properties of the Ising model. In fact, a number of efforts have been 
preformed on the Ising model with lattices whose topology is equivalent to a spherical 
surface |in | ITT ] IT2 \ I18 | IT^ . It is noted that, however, the thermodynamic limit cannot 
be considered for the closed form of sphere-like surfaces having positive curvature while 
maintaining their finite curvature. This is because a spherical surface reduces to a flat 
plane in this limit, where the curvature effect vanishes completely. 

Therefore, instead of a sphere, we consider a curved surface with negative constant 
curvature, termed a pseudosphere [THl EDj- The pseudosphere is a simply connected 
infinite surface in which the Gaussian curvature at arbitrary points possesses a constant 
negative value. (The definition of the pseudosphere will be given in Appendix.) Hence, 
it serves as a suitable geometry for considering the curvature effect on the critical 
properties of a system. It should be noted that the pseudosphere occurs in manifold 
physical problems ranging from quantum Hall effects [211 EHl IIEI 1211 12H| , quantum chaos 



Geometric effects on critical behaviours of the Ising model 



5 



X 10' 
8 I ^ 




Linear dimension of lattices r 



Figure 2. Total number of sites N involved in the hcptagonal lattice with various 
system sizes. The horizontal axis represents the number of concentric layers of 
heptagons, r; this effectively serves as the linear dimension of the entire lattice. Inset: 
A single logarithm plot of the identical data shown in Fig. 2. It is clearly seen that N 
for r 3> 1 increases exponentially with r. 

j2niEniEIl5 the string theory [22] to cosmology [221, wherein the underlying geometric 
character of the system is extremely significant. 

Interestingly, a wide range of regular lattices can be constructed on the 
pseudosphere [201 • This is achieved by a tessellation procedure, where the entire surface 
is covered by non-overlapping regular polygons meeting only along complete edges or 
at vertices. It is known that a regular tessellation of the pseudosphere with q regular 
p-sided polygons meeting at each vertex satisfies the following property: [20] 

(p-2)(g-2)>4. (4) 

Hence, the series of integer sets {p, q} satisfying (jlj results in an infinite number of 
possible regular tessellations of a pseudosphere. This is in contrast to the case of a 
flat plane, where only three regular tessellations are allowed: {p, q} = {3, 6}, {4, 4} and 
{6, 3} satisfying the condition (p — 2)(g — 2) = 4. For simplicity, we adopted a heptagonal 
{7, 3} tessellation to construct the Ising lattice on a pseudosphere. Figure 1 illustrates 
the local bond structure of a regular heptagonal lattice in terms of the Poincare disk 
representation. The resulting lattice comprises concentric layers of congruent heptagons 
surrounding a central heptagon. The Ising lattice models embedded on a pseudosphere 
have been considered thus far [211 EH ESI ; however, explicit temperature dependences 
of thermodynamics quantities close to the transition are yet to be concerned. 

Due to a peculiar metric of a pseudosphere, the total number of sites of our 
heptagonal lattice exhibit a non-trivial evolution behaviour with the increase in the 
lattice size. The size of our lattice is determined by the number of concentric layers of 
heptagons, denoted by r, which effectively serves as a linear dimension in our lattice. 
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For a given r, the number of total sites is expressed as follows: 
N{r = 1) = 7, 

r-2 

N{r > 2) = 7 + 7^^ 

where c± = 2±v^- Figure|21plots the dependence of on the effective linear dimension 
r. When r ^ 1, it is approximated as N{r) ~ 5exp(r); this means that N rapidly 
increases with r in comparison with the case of the planar Ising model. The exponential 
increase in N{r) is a manifestation of the constant negative curvature of the underlying 
geometry of our lattice. 

It should be noted that, when considering thermodynamic properties of our lattice, 
careful treatments on boundary effects are required. The exponential increase in N{r) 
results in that the ratio [N{r) — N{r — 1)]/N{r) approaches a non-zero constant 1 — e^^ 
in the limit r ^ oo. This means that the boundaries of our lattices can not be 
neglected even in the thermodynamic limit, but contain a finite fraction of the total 
sites. Boundary effects coming from these sites are difficult to be eliminated completely, 
because the periodic boundary conditions are hard to be employed to regular lattices 
assigned on a pseudosphere. 

In order to extract the bulk critical phenomena, therefore, we have followed the 
procedure mentioned below. Suppose that an Ising lattice consists of rout concentric 
layers of heptagons. Then, for computing physical quantities of the system (magnetic 
susceptibilities, for instance), we take into account only the Ising spins involved in the 
interior Tin layers (rin < rout) so as to reduce the contribution of the spins locating near 
the boundary. In actual calculations, Tin is varied from 4 to 8, and for each Tin the 
number of disregarded layers Ar = rout — ''"in is systematically increased from to 4. 
By investigating the asymptotic behavior of the system for large Ar, we can deduce the 
bulk properties of the Ising lattice model embedded on the pseudosphere. 



•1 + C+V 



l + c_ 



(5) 



4. Numerical methods 



We considered the conventional Ising model with ferromagnetic interaction: 

H = -J SiSj, Si = ±1, (6) 

<i,j> 

where denotes a pair of nearest-neighbour sites on a heptagonal lattice. The free 
boundary condition is imposed for all lattices to be considered. Temperatures and 
energies are expressed as units of J/ /cb and J, respectively. The order parameter m per 
site for a given configuration of {sj} is given by m = Y^iLi Si/N. 

Our main objective is to determine the zero-field magnetic susceptibility x that 
exhibits the power-law relation x{T) oc |T — Td""^ near the critical temperature T^. 
The susceptibility for a finite system size can be expressed as a function of the order 
parameter m as 

X(T,.,.^, (7) 
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or alternatively, 



x'(T,iV) = iV 



(m^) - {\m\f 



(8) 



Despite the difference in their definitions, both x and x' yield the same critical exponent 
7 by the finite-size scaling method, as described in Ref jHEl- The expectation values 
(|m|) and (m^) at temperature T are evaluated by canonical MC simulations [5^1 137j. 
Sampling of the configurational space was carried out by using the single-cluster update 
algorithm which prevents the critical slowing down near the transition. 

Quantitative evaluation of critical exponents can be achieved by using the finite-size 
scaling technique [SHI EDI E] ■ Close to the critical temperature Tc, the susceptibility for 
a finite system size satisfies the following scaling behaviour: 



Here, /i is the critical exponent describing the divergence of the correlation volume ^v(T) 
of the order parameter: 



The quantity is a natural generalization [23 112] of the correlation length ^ that 
diverges as C,{T) oc \T — Tc\~'^ with the critical exponent u in the planar model. Near T^, 
the argument of the scaling functions Xo in ©, denoted by x = |T — Tc\N^^^, becomes 
much smaller than unity. This allows the polynomial expansion of the scaling function 

Xo as 



terminating the expansion at the order n. By substituting the numerical data of 
x{T, N) and their corresponding values of T and into (jllll . followed by performing the 
nonlinear least square fitting, we obtain the critical exponents 7 and /i and the critical 
temperature as optimal fitting parameters. 

5. Results 

5.1. Susceptibilities and critical exponents for Ar = 

Before addressing the bulk critical properties, we first demonstrate the results for entire 
heptagonal lattices with Ar = (boundary contributions are fully involved). Figure 
El^a) andin^b) show the calculated results of the zero-field susceptibilities x{T,N) and 
x'{T,N), respectively, as a function of temperature T. The single logarithmic plot is 
used in the figure IHl^a). The system size r = Tin is varied from 4 to 8, which corresponds 
to the change in the number of total sites from = 315 to = 15435. Both x and 
x' exhibit such a typical behaviour that indicates the occurrence of a ferromagnetic 
transition within the temperature range of 1.1 < T < 1.3. For instance, the x curve in 
figure infa) monotonically increases with a decrease in T; this is attributed to the onset 




(9) 



ev(T) cx IT - Te 



(10) 



X{T,N) = aoN^'^ + ai\T- T^\N^'^+^^/^' + ■■■ + a„|T - Tc|Ar("+^)/^, (11) 
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Figure 3. (a) Zero-field magnetic susceptibilities: (a) x{T,N) and (b) x'{T,N), for 
entire heptagonal lattices with Ar — 0. The system size is varied from r-m — 4 (square) 
to Tin = 8 (solid circle), which corresponds to the change in the number of total sites 
from TV = 315 to TV = 15435. 

of the ordered phase. As well, the plot of x' in figure E^b) exhibits a sharp peak at 
T ~ 1.2, which is also a precursor of the divergence in the infinite system. 

Figure El (a) shows the scaling plot of x(T,N) based on Q. The vertical 
and horizontal axes represent the scaled susceptibility xi^)^^'^^'''' ^ind its argument 
X = \T — Tc\N^^^, respectively. The critical exponents evaluated from the upper and 
lower branches are 7 = 2.28(2) and /i = 3.46(1), and 7 = 2.26(2) and /i = 3.47(1), 
respectively. The errors in the last decimal places, which are shown in parentheses, 
designate a 95% confidence interval. As expected, the estimated values of 7 and /i for 
the two branches are in agreement within numerical errors. The optimal value of the 
critical temperature for the two branches is evaluated as Tc = 1.253(1); this agrees with 
the preceding estimation from the values of figures a) andEfb). A similar analysis 
for x' results in 7 = 2.26(3), /x = 3.45(2), and = 1.254(2); these values are fully 
consistent with the results deduced from the data of x- 

We re-emphasize the fact that for regular Ising lattices embedded on curved 
surfaces, instead of r should be adopted as the scaling variable. This is justified 
by attempting the finite size scaling of x using r. Figure EJ^b) presents the scaling plot 
of X terms of another scaling argument y = \T — Tc\r^^^'; the estimated values 
of the parameters are 7' = 1.8(3), fi' = 2.0(4), and Tc = 1.46(4). It is evident 
that the data points do not collapse onto a single curve, but instead exhibit a large 
scatter. Furthermore, the resulting value of Tc differs from the temperature at which 
the susceptibility x'(T) exhibits a sharp peak (see figure 01). These facts indicate that 
the linear dimension r is not a characteristic length scale that describes the scaling 
behaviour of the thermodynamic quantities for curved surfaces. 
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Figure 4. (a) Scaling plot of xiT,N) against the argument x = \T — Tc\N^^f^. 
The optimal values of critical exponents are 7 = 2.28(2) and /i ~ 3.46(1), and 
7 — 2.26(2) and /x = 3.47(1) for the upper and lower branches, respectively. The 
critical temperature is estimated as Tc — 1.253(1) for the two branches, (b) Scaling 
attempts of x against the alternative argument y = jT — Tc|r^/'' instead of x. Optimal 
values of parameters are 7' = 1.8(3), /i' = 2.0(4), and Tc = 1.46(4). 



5.2. Critical exponents for Ar > 0; Boundary effects 

We now turn to the study of bulk critical properties of the heptagonal Ising lattice 
model. As mentioned in Sec. 3, the boundary spins of Ising lattices on a pseudosphere 
are thought to affect significantly to the nature of the system, since the number of the 
spins along the boundary increases as fast as that of total spins of the lattice. Hence, 
in order to extract the bulk critical exponents, we must try to remove the contribution 
of the boundary spins to the scaling behavior of the system. This is achieved by setting 
the disregarded layers Ar = rout — ''"in to be finite, i.e., by summing up only the spins 
within the interior r[^ layers when performing MC simulations on the systems consisting 
of '"out(> '"in) layers. If Ar is sufficiently large, the ensemble of the spins involved in the 
interior rjn layers may yield the critical exponents of the bulk system that is free from 
the boundary contribution. 

On the basis of the argument above, we have prepared the heptagonal lattices 
having various values of r^^ and Ar, and systematically employed the scaling analysis to 
them. It then reveals how the critical exponents 7 and /i and the critical temperature 
Tc depend on the number of disregarded layers Ar. The calculated results are given in 
Figure El (a) and El (b); each data point in the plots was extracted by means of the finite 
size scaling analysis for the system sizes 4 < r^n < 8. (Hence, the maximum system 
size we have treated reaches rout = 12 which corresponds to iV = 725760.) We found 
that an increase in Ar results in a monotonic decrease in all quantities in question: 7, /i 
and Tc, which indicates the significant contribution of the boundary spins to the scaling 
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Figure 5. (a) Ar-dependences of the critical exponents 7 and /i, and (b) tliat of tlie 
critical temperature T^. The lines in the plots serve as a guide to eye. Each data 
point was extracted by means of the finite size scaling analysis for the system sizes 
4 < Tin < 8. The values of the mean-field exponents 7mf and /imf are indicated in the 
plot (a). 



behavior of the system. 

Asymptotic behavior of the curves of 7 and fi for large Ar provide estimations 
of the bulk critical exponents. From Figure |S1 (a), we see that the curve of /i shows 
somewhat convergence to a particular value of ~ 2 or less. Since this is close to the 
value of /i for the two-dimensional planar Ising model||, fi2d = 2, it appears that the 
bulk system attains the planar Ising universality class at around Ar ~ 4. However, the 
asymptotic behavior of 7 is slightly different from that of /i; while the value of 7 equals 
to that of the planar system 72d = 1-75 at Ar ~ 3, it still continues to decrease almost 
linearly with Ar and thus has no tendency to converge to 72d for large Ar. Thereby, the 
bulk critical exponent 7 will be smaller than 72d, which suggests that the heptagonal 
Ising lattice belongs to some other universality class than of the planar Ising lattices. 

In the context above, the asymptotic value /i ~ 2 for large Ar is not the exponent 
fi2d but another specific exponent characterizing the intrinsic nature of the system. 
This point is clarified by referring to the previous studies done by Rietman et al. j2I] 
and by Doyon and Fonseca |22]- They both have stated that the Ising lattice model 
embedded on the pseudosphere should yield the mean-field critical exponents when 
the boundary contribution may be omitted. The mean-field nature of the system is 
attributed to the fact that an Ising lattice embedded on a pseudosphere is effectively 
an infinite dimensional lattice at large distance due to the exponential growth of the 
total spins 132] ■ For an ordinary Ising lattice in d dimension, the number of spins along 
the boundary, iV^, is related to that of the total spins as A^, oc A^i-(Vrf). Hence, 

11 The value ^2d = 2 is derived from the relation ji — vd (See and the exact solution = 1 for 

d^2. 
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the peculiar relation Ns (x N that holds on negatively curved surfaces consequences 
d = oo. Accordingly, our heptagonal Ising lattices are expected to yield the mean-field 
critical exponents 7mf = 1 and /imf = 2, where /zmf is determined by^ /imp = T-'Mpdc as 
suggested in Ref. |23|. Our numerical results given in Figure El (a) is in fact qualitatively 
consistent with the argument above; that is, converges to /imf = 2 for Ar > 4, and 
7 continues to decrease until it yields 7mf = 1- (To be precise, the asymptotic value of 
7 may take the value between 72d and 7mf depending on the numerical conditions; this 
point will discussed in Section 6 in detail.) 



5.3. Binder's cumulant U4{T, N) 

We have also carried out an alternative estimation of Tc and n in terms of the fourth- 
order Binder's cumulant Ui(T,N) defined by [321 EH Ell 

U.iT.N) = l-^^. (12) 

For a given A^, the cumulant Ui{T) decreases monotonically with an increase in T from 
Uiiff) = 2/3 to U^^oo) = 0. In the vicinity of the T-dependence of it can be 
approximated by 

f/4(r, A^) = [/W + [/(I) (^1 - AtV/^, (13) 

where U^^^ and U^^^ are constants and are thus independent of T and A^. The expression 
()13|) is a direct consequence of the assumption that the probability distribution of the 
order parameter m should be Gaussian close to the transition [3Q\. From ()13|1 . it follows 

dU4 



UiiTc.N) = const, and 



dT 



oc -A^^/^, (14) 



which provides a complementary method to estimate Tc and /i. 

Figures El (a)-(c) present the numerical results for Ui{T,N) for several values of 
Ar. In each plot, we found a unique crossing point giving an estimate of the critical 
temperature as Tc ~ 1.25, 1.22 and 1.20 for Ar = 0, 2 and 4, respectively. These 
values of are in fair agreement with those obtained by the scaling analyses for the 
susceptibilities presented in Figure El (b). 

FigureEl(d) shows the A^ dependence of the derivative with negative sign —dU^^/dT 
at T = Tc for different Ar. The magnitude of the derivative grows with a power-law 
with increasing A^ as expected from (|T^ . The estimates of /i are /i ~ 3.6, 7.7, 12.5 
for Ar = 0, 2, 4, respectively. It should be noted that these values of fi increase with 
Ar. This clearly contradicts to the results of the scaling analyses on the susceptibility x 
(See Fig. El (a)), where /i is found to be such a decreasing function of Ar that yields the 
mean- field value /imf = 2 in the large Ar limit. The discrepancy may be due to simply 
a finite sized effect; or, it may indicate some intrinsic property of negatively curved 

% For the Ising model, the mean-field exponent vmf is 1/2, and the upper critical dimension dc is 4; 
thus we obtain /imf = 2. 
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Figure 6. (a)-(c) Temperature dependences of the fourth-order Binder's cumulant 
U^lTjN) for various values of Ar. The size of the interior layers is varied from 
Tin — 4 (square) to rin = 8 (solid circle), (d) The negative slope of the cumulant, 
—dUi/dT\T=Tc, for each value of Ar. The solid lines represent the power-law (c>c iV^/'^) 
with fi — 3.6, 7.7, and 12.5 from top to bottom. 



surfaces with regard to the distribution of the order parameter m, since the power-law 
relation ()14|1 originate from the assumption of the Gaussian distribution of m j2E] • More 
through discussion about this point will be presented elsewhere. 

6. Discussions and Concluding remarks 

Our numerical analysis revealed that the critical exponents 7 and fi of the heptagonal 
Ising model defined on negatively curved surfaces assume values that deviate from those 
for the planar Ising model. Most striking is that, when reducing the contribution of spins 
near the boundary, both 7 and /i exhibit a tendency to yield the mean-field exponents. 
This phenomenon is attributed to the fact that the regular Ising lattices on negatively 
curved surfaces serve as an effectively infinite dimensional lattices due to the peculiarity 
of the intrinsic geometry. 

The above statement immediately poses the following question: Does the negative 
curvature of the underlying geometry alter the other four critical exponents? With 
regard to the power-law behaviour of thermodynamic quantities, the planar Ising model 
is known to possess four other critical exponents a = 0, /3 = 1/8, 5 = 15, 
and rj = 1/4, which correspond to heat capacity, spontaneous magnetization, critical 
isotherm, and the two-point correlation function, respectively. Our preliminary study 
jlH] has suggested that, the exponent (3 for large Ar also tend toward the mean-field 
exponent (3mf = 1/2, whereas Tc estimated there is slightly different from that in the 
present work. Detailed analyses on this issue and the quantitative determination of the 
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other exponent 5 (and 77, if it exists^) will be given in a future study. Very recently, 
we have found that the dynamic critical exponent z in our heptagonal Ising model also 
shift quantitatively from that for the planar Ising model ^H] . 

It deserves comment that the finite curvature of the underlying geometry may 
produce another type of effect wherein the spin variables at each site possess 
orientational degrees of freedom. This is because the relative angle of interacting spins 
at neighbouring sites on a curved surface is determined by a spatially-dependent metric 
tensor. Thus, the Hamiltonian of the system should be modified such that it is a 
function of the metric tensor. As a result, the energetically preferable configurations of 
the vector spins differ from those in planar systems jTft 0HI UHl EDI which implies 
that the critical behaviour of these vector-spin lattice models is markedly influenced by 
a finite surface curvature. 

Further noteworthy is, however, the fact that the geometric curvature continues 
to be relevant to the critical behavior of the system despite the omission of the 
vector property of the interacting spin variables; this was demonstrated by our results. 
Evidently, the spin variable of our model is set to be a scalar, and thus the Ising 
Hamiltonian given by (jH)) is devoid of the metric factor. Nevertheless, the surface 
curvature is surely relevant to the Ising model on curved surfaces since it enables to 
construct peculiar lattice structures that cannot be realize in a flat plane. This means 
that the surface curvature induces the alteration of the global symmetry of the system 
even when the interacting entities do not exhibit any vector property. Besides, in the 
vicinity of the critical point, the discreteness of the lattice becomes irrelevant and the 
model can be considered to be a continuum surface with a constant curvature. Therefore, 
the metric of the underlying geometry plays a crucial role in the scaling behaviour of 
the Ising lattice model defined on the surface. In this context, the mean-field nature 
of the critical exponents 7 and /i are expected to be universal for all lattice structures 
other than the heptagonal one; this point is being investigated. 

We remark that it is also interesting to study the dependence of the values of the 
bulk critical exponents on the interior lattice size Tin we have introduced. Obviously, 
the condition 1 -C <C rout is desirable to determine the bulk properties of the lattice 
with accuracy. If the rin is not so large (compared to the curvature radius of the 
underlying surface), the system can be regarded as an Ising lattice defined on a nearly 
fiat surface. Thereby, the resultant critical exponents will become comparable with 
those for the planar Ising lattice rather than the true bulk critical exponents (i.e., the 
mean-field exponents). This implies that, by gradually increasing the size of rjn, the 
system go through a crossover from the planar Ising class to the mean-field class. We 
conjecture that this crossover phenomenon is observed in the shift of the asymptotic 
value of 7 deduced from the plot in Figure El (a); that is, when increasing the size of 
r employed in the analyses, the asymptotic value of 7 at Ar ^ 1 will shift downward 

+ Rietman et al. [23 have suggested that the two-point correlation function of the Ising lattices on 
a pseudosphere shows an exponential decay instead of an power-law decay. If it is true, the critical 
exponent 77 can be no longer defined. 
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from 7 ~ = 1.75 to 7 ~ 7mf = 1- Accordingly, it is possible that the curve of 7 
plotted in Figure El (a) converges to an intermediate value between 72d and 7mf under 
the present numerical conditions. 

In conclusion, we have investigated the critical behaviour of the Ising model defined 
on a curved surface with constant negative curvature. MC simulations and finite-size 
scaling analyses were employed to compute the critical exponent 7 and n for the zero- 
field magnetic susceptibility and correlation volume, respectively. The resulting values 
7 and /i show distinct values from those for the planar Ising model, and exhibit a 
tendency to the mean-field exponents 7mf = 1 and /imf = 2 due to the peculiar intrinsic 
geometry of the negatively curved surface. As well, we have revealed quantitatively how 
the boundary spins contribute to the determination of 7, fi and T^, and argued the 
possibility to occur the crossover from the planar Ising class to the mean-field Ising 
class. We hope that the generalization our statistical model (with regard to the lattice 
structure, dimensionality of the embedding space, distribution of interacting strength, 
etc.) would unveil a wide variety of interacting critical properties of the physical systems 
assigned on general curved spaces. 
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Appendix A. The pseudosphere 

The pseudosphere is defined as one sheet of the double-sheeted hyperboloid 

x'^ + - = -1, (A.l) 

possessing the Minkowskian metric ds"^ = dx^ + dy^ — dz^. Since (jA.lj) specifies the 
locus of points whose squared distance from the origin are equal to —1, it is called a 
pseudospheie having the radius i = yj—l by analogy with the sphere. 

While the above definition is rigorous, it is clumsy for computations since three 
coordinates are used for only two degrees of freedom. This cumbrousness is removed by 
using an alternative representation of the pseudosphere, called the Poincare disk model 
j2n]- Suppose that the upper hyperboloid sheet is projected onto the x-y plane using 
the following mapping: 

f-^, (A.2) 



,1 + 2 l + z 

This transforms the upper sheet to a unit circle endowed with the metric 

2 

ds^ = w"^ idx^ + dy'^) , w = -. (A. 3) 

1 — x'^ — 



The unit circle possessing the metric flA.3|) is referred to as a Poincare disk, and it serves 
as a compact representation of the pseudosphere. The boundary of the disk corresponds 
to the points at infinity of the hyperboloid. The Gaussian curvature k on the disk is 
calculated using the formula: 

1 / 92 52 \ ^ ^ 

w"^ \dx^ dy"^ ) 

From ()A.3|) and ()A.4|) . we see that k = — 1 at arbitrary points on the disk. Thus it 
follows that the pseudosphere is a curved surface with a constant negative curvature 
K = -1. 
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